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We present a computation of the 0{atas + a?) two-loop corrections to the MSSM Higgs 
masses. An appropriate use of the effective potential approach allows us to obtain simple 
analytical formulae, valid for arbitrary values of niA and of the mass parameters in the stop 
sector. In a large region of the parameter space the 0{a^) corrections are comparable to the 
0{atas) ones, increasing the prediction for mn by several GeV. 



In the Minimal Supersymmetric Standard Model (MSSM), the Higgs sector contains two 
SU{2)l doublets, and the masses and the couplings of the Higgs bosons can be expressed, at the 
classical level, in terms of known Standard Model quantities and only two additional parameters: 
for example, the mass of the neutral CP-odd state, m^, and the ratio of the two Higgs vacuum 
expectation values, tan j3. In particular, the mass rrih of the lightest Higgs boson is predicted to 
be smaller than the mass of the Z boson. However, radiative corrections modify considerably 
the predictions for the Higgs masses, allowing the light Higgs mass to escape the tree-level 
bound, and bring along additional dependences on the remaining MSSM parameters. Direct 
searches at LEP have already ruled out a considerable fraction of the MSSM parameter space, 
and the forthcoming high-energy experiments at the Tevatron and the LHC will either discover 
(at least) one light Higgs boson or rule out the MSSM as a viable theory for physics at the weak 
scale. Thus, a precise calculation of the radiative corrections to the Higgs boson masses, valid 
for arbitrary values of all the relevant parameters, is needed to compare reliably the MSSM 
predictions with the (present and future) experimental data. 

In this talk, a computationlllfl of both 0{atas) and 0{a'i) two-loop corrections to the MSSM 
Higgs boson masses, based on the effective nptential approach, is presented. The results of refs.EIS 
extend the existing two-loop calculations S~Q to the case of arbitrary MSSM parameters, and 
are given as a set of simple analytical formulae, suitable for being promptly included in the 
experimental analyses. 

We start our discussion of the corrections to the Higgs masses by recalling that in the 
effective potential approach the CP-odd and CP-even mass matrices are identified with the 
second derivatives of the effective potential, 14ff = Vq + Al/, evaluated at its minimum: 



where we have decomposed the Higgs fields into their VEVs plus their CP-even and CP-odd 
fluctuations as = {vi + Si + i Pi) / \/2 . Using the explicit expression of the tree-level potential 
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Vo, and imposing minimization conditions on Ves, the mass matrix for the CP-even sector can 
be written as: 
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Ml) = {Ml) + {ami) > (2) 

where 

y - {ml + ml) sp cp m| + j ' ^ ' 
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In Eq. (|^, ml is the mass of the CP-odd state in the effective potential approximation, i.e. 
the non-vanishing eigenvalue of [M^pY^ . Similarly, m| = {g^ + g'^)f^/4 is just the DR mass 
for the Z boson. Finally, sp = sin (3 and c/3 = cos /?, where tan (3 = V2/V1. 

According to Eq. @, (AA^I)*^^ can be computed by taking the derivatives of AV with 
respect to the CP-even and CP-odd fields, evaluated at the minimum of 14fj. In this computation 
we express the field-dependent masses and interaction vertices that contribute to Ay, and are 
relevant for our calculation, in terms of five field-dependent quantities. The first three, appearing 
already in the one-loop contribution to the effective potential, are the squared masses of the top 
quark and of the stop squarks. The remaining two quantities, appearing only in the two- loop 
contribution to the effective potential, can be chosen as the stop mixing angle, denoted as and 
cos {(p — (p), where if is the phase in the complex top mass, and ip is the phase in the off-diagonal 
term of the stop mass matrix. 

After a straightforward although lengthy application of the chain rule for the derivatives 
of the effective potential, we derive an expression for the two-loop corrections to the Higgs 
mass matrix, valid when the one-loop part of the corrections is written in terms of one-loop 
renormalized parameters: 

{ami) = "'^ ^/Z , (5) 

where Z is a wave function renormalization matrix for the Higgs fields, and 

\ efF 1 

AMI) = 2 + 2 A.Fs) , (6) 

AMiy^^ = hi /i mt S26 (F2 + A^Fs) + ^ At ^ sle (F3 + Aa.Fj + A^F^) , (7) 

AMiy^^ = 2 h?t ml Fi + 2 hi At mt S2e {F2 + A^.Fs) + ^ hi Al sle (^3 + 2 Aa,F^) . (8) 

In the above equations the functions (^1,^2,^3) are written as Fi = Fi + AFj, where Fj are 
combinations of the derivatives of the pure two-loop effective potential whose explicit definitions 
are given in ref.§. The terms AFj include contributions coming from the renormalization of the 
top and stop masses appearing in the one-loop part of the Fj, as well as from the renormalization 
of the common factors multipliying Fj, i.e. {hi ,mt ,S2e)- On the other hand, in each entry of 

{AMiy the functions Fj are multiplied by different combinations of // and At. These two 
parameters have different 0{at) and 0{as) renormalizations that cannot be absorbed in the F,, 
and are then separately taken into account by the terms A^Fj and A^i^Fj. 

We have derived explicit analytical formulae for the 0{atas) and 0{al) parts of the functions 
Fj, obtained under the assumption that the 0{at) parts are expressed in terms of DR quantities. 



In ref. we presented the complete formulae for the 0{atas) corrections, valid for arbitrary 
choices of all the input parameters. On the other hand, the general formulae for the 0{a1) 
corrections are quite long, thus we made them available upon request in the form of a computer 
code. In ref. we presented explicit analytical formulae valid for mq = mu = AIs, in which 
0{mt/Ms) corrections are neglected, but is left arbitrary. 

In the phenomenological analyses of the MSSM Higgs sector, it may be useful to express 
or results in terms of input parameters given in a different renormalization scheme, which will 
be indicated generically as R. To obtain the 0{atas + af) correction in the R scheme one has 
just to shift the parameters appearing in the one-loop term, i.e. = + 6x, where x 

is a generic parameter. This induces additional two-loop contributions that can be absorbed 
in a redefinition of the functions Fi and AFi entering Eqs. Explicit formulae for the 

transition from the DR scheme to a generic R scheme are given in refs. ElS. 

The parameters entering the one-loop part of the corrections can be identified as {rrit, , 
m^-^, At, fj,, V, tan/3). Some of them are related by the equality 

sin2%= '"^^^t + ^f (9) 
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We choose to express our result identifying the stop masses rrij and the top mass rrit with the 
corresponding pole masses. For the electroweak symmetry-breaking parameter v, we use the 
value obtained in terms of the precisely known muon decay constant G^. These choices specify 
6mf , 6mt and 5v in terms of the finite parts of stop, top and W-boson self-energies respectively. 

ti 

For the stop mixing angle, several OS prescriptions are present in the literature. We find suitable 
for our 0{atas + a1) calculation the following 'symmetrical' definition^ 

1 ni2(m? ) + ni2(m? ) 
Sdr = -Ay , (10) 

where ni2(p^) stands for the finite part of the off-diagonal stop self-energy. 

It is not clear what meaning should be assigned to an OS definition of the remaining quan- 
tities, i.e. (Aj, tan /3). For instance, they could be related to specific physical amplitudes. 
However, given our present ignorance of any super symmetric effect, such a choice does not seem 
particularly useful. Then, we find simpler to assign the quantities (i and tan (3 in the DR scheme, 
at a reference scale Qq = 175 GeV, that we choose near the present central value of the top 
quark mass. In this framework, we are going to take 6^ = 6[3 = 0, and we treat At as a derived 
quantity, obtained through Eq. 

Our effective potential calculation is equivalent to the evaluation of the Higgs self-energies 
in the limit of vanishing external momentum. A diagrammatic computation of the two-loop 
0{atas) contributions to the Higgs boson self-energies at zero external momentum has been 
performed in ref. §. Analytical formulae, valid in the simplified case of degenerate soft stc 
masses and zero mixing (with /i = At = 0), have been presented in the first paper of ref. 
For arbitrary values of the top and stop parameters, however, the complete analytical result of 
ref. □ is far too long to be explicitly presented, and is only available as a computer codeB. We 
have checked that^ in the case of zero mixing and degenerate stop masses, our results coincide 
with those of ref. Moreover, after taking into account the difference in the definitions of the 
OS renormalized angle 9i, we find perfect agreement with the numerical results of ref. for 
arbitrary values of all the input parameters. 

A calculation of both 0{atas) and 0{a'^) two-loop corrections to the lightest Higgs boson 
mass m/i, based on the formalism of the effective potential, has been presented in refs. In 
these papers, however, the dependence of the stop masses and mixing angles on the fields Pi 
(the CP-odd components of the neutral Higgs fields) is not taken into consideration. Therefore, 
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Figure 1: The mass rrih as a function of ttia, for tan (3 = 2 or 20 and Xt = or 2 TeV. The other parameters 



are rrin = rn. 



= 1 TeV, n = 200 GeV, nig — 800 GeV. The meaning of the curves is explained in the text. 



the OioitOLs + a^) corrections to the parameter mA are not evaluated. If one wants to relate 
the input parameters to measurable quantities, the computation is applicable only in the limit 
rtiA ^ Jn-Z) in which nih is nearly independent of m^. Moreover, the results of refs. 10 for rah 
are available in numerical form, and simple analytical formulae are provided only in the limit of 
universal soft stop masses (degenerate with txia and m^) much larger than the top quark mass. 
We have verified that, in such limit, our analytical results for rrih agree with those of refs. 

In the analysis of ref. § some 'leading logarithmic' 0{a^^ corrections, obtained by renormal- 
ization group methods have indeed been added to {M'^s)22- However, when comparing our 
complete O(a^) result with the renormalization group result, it turns out that, for some choices 
of the SUSY parameters, the leading logarithmic corrections amount only to a fraction of the full 
ones. Fig. || shows m/j as a function of ttia, for tan /? = 2 or 20 and|| X'^'^ = or 2 TeV. The other 
input parameters are chosen as rat = 175 GeV, raQ^ = = 1 TeV, /i = 200 GeV, mg = 800 
GeV. The cases with X^^ = and X^^ = 2 TeV correspond, respectively, to the so-called 
'no-mixing' and 'm/j-max' benchmark scenarios t3 considered in the experimental analyses EJ. 
The curves in Fig. |^ correspond to the two-loop corrected Higgs mass at 0{atas) (short-dashed 
line), at 0{atas + af) including only the leading logarithmic corrections (dot-dashed line), and 
at 0{atas + a^) including our full computation (solid line). The one-loop result for rah is also 
shown for comparison (long-dashed line). It can be seen from the figure that the 0{atas) cor- 
rections are in general large, reducing the one-loop result for mh by 10-20 GeV. On the other 
hand, the 0{af) corrections tend to increase nih- for small stop mixing they are generally small 
(less than 2-3 GeV), whereas for large stop mixing they can reach 7-8 GeV, i.e. a non-negligible 
fraction of the 0{atas) ones. Moreover, it appears that the 'non-leading' corrections included 
in our full result are always comparable in size with the 'leading logarithmic' ones, and are even 
more significant than the latter in the case of large stop mixing. For ^ ruz, we find good 
numerical agreement with the value of m/j obtained from the approximate formulae of ref. Q. 



"To facilitate the comparison with the existing analyses of the two-loop corrected Higgs masses, we make use 
of the unphysical parameters (mQ^,m^^, X?^) that can be derived by rotating the diagonal matrix of the OS 
stop masses by an angle Sj. 



In conclusion, refs.u'B contain the most advanced calculation (so far) of the two-loop cor- 
rections to the MSSM neutral Higgs boson masses. Our work should lead to a more accurate 
interpretation of the experimental searches for the neutral MSSM Higgs bosons at LEP, the Teva- 
tron, the LHC and other possible future colliders. The importance of the new 0{a^) effects we 
have computed will increase further when the top quark mass is measured more precisely. Then, 
we can hope for the next step, the discovery of super symmetric particles and super symmetric 
Higgs bosons. 
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